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1. INTRODUCTION 

The primary objective of Energy control centre is to maintain the power system in a secure and 
stable state by continuously monitoring the power flows in the lines and voltage magnitudes at the buses. 
Voltage variations are due to the imbalance of reactive power generated and consumed by the node. These 
voltage variations can be corrected by co-ordinated control of voltage/reactive power control devices such as 
generator excitations, on-load tap changing transformers and switchable shunt VAR compensating devices. 
The complexity of power system operation is increasing day by day because of growing demand and without 
matching generation and transmission facilities resulting transmission as well as generation operated with 
smaller safety margins. Under these conditions, it is great challenge to optimize the power system and ensure 
the security. 

Optimal Reactive Power Dispatch (ORPD) problem is a non linear optimization problem with 
multiple objectives and constraints. This problem has been solved by a number of conventional optimization 
techniques such as Linear programming (LP), Non linear programming (NLP), Quadratic programming etc. 
are reported in the literature [1-6]. However, these conventional algorithms have some limitations such as 
differentiation of objective function is required and curse of dimentionality. These limitations can be 
overcome if evolutionary computation techniques are adapted because of their approach of random search 
and begin with a population of solutions and also no differential information is required. There are various 
evolutionary computation techniques such as Particle Swarm optimization algorithm, Gravitational search 
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algorithm, Firefly algorithm, Differential evolution algorithm, ant colony algorithm, BAT algorithm etc. are 
reported in the literature [7-11] for ORPD problem. 

Bacterial Foraging Algorithm (BFA) is one of the population based evolutionary computation 
algorithms which is proposed by Passino. This algorithm is inspired by foraging strategy of E.coli bacteria. 
Though original BFA algorithm is used by many experts in different fields, it has some drawbacks of slow 
convergence speed and not taken into account on the diversity of population which can easily lead to the 
premature convergence when dimension increases. Yuan —-tao Zhang et al. [12-13] proposed adaptive 
chemotaxis step setting and chaotic perturbation in each chemo tactic to improve the original BFA. Fuqing 
Zhao et al.[ 14-15] proposed a chmotaxis enhanced BFA by introducing differential mutation operator and 
chaotic operator in chemotaxis step. Na Dong et al.[16] proposed chaotic PSO algorithm by introducing 
chaos dynamics to generate initial population and chaotic perturbation in to swarm updation. 

The present paper proposes a Modified Bacterial Foraging algorithm (MBFA) which utilizes the 
differential mutation operator to enhance the tumble step in chemotaxis to overcome tumble failure and chaos 
dynamics to generate uniformly distributed initial population. The proposed algorithm is applied to solve the 
optimal reactive power dispatch with two objectives minimization of real power loss and voltage stability L- 
index. The performance of the proposed algorithm is tested on IEEE 30 bus standard test system and results 
are compared with other evolutionary techniques available in the literature on this problem. The superiority 
of proposed algorithm is demonstrated. 


2. PROBLEM FORMULATION 
The proposed algorithm is applied to two objective functions: minimization of real power loss and 
voltage stability L-index. 


2.1. Transmission loss objective (Ploss): 
The real power loss of the system can be calculated as follows 


Nline 
Poss = Bet ox Ve a VP — 2ViV; cos(6; mi 9;)) MW (1) 


Where Pioss is the total real power loss, Njine is total number of transmission lines. Vjand V; are the voltage 
magnitudes at the two ends of the K" line. 6; and 6; are the voltage angles at the two ends of the K” line. 
Voltage magnitudes and angles can be caluculated from the load flow solution. g; is the conductance of the 
K"" line. 


2.2. Voltage stability index objective (Vstability): 
Voltage stability L-index is considered as measure to find voltage stability. 


Linax= max(L)) (2) 


L-index can be computed as 
L=l-) FF. (3) 
V. 


Where j indicates all the load buses. v; and v; are voltage magnitudes at i" and j" buses and they can be taken 
from load flow. Fji can be obtained from the Ybus matrix as follows 


I Yic Yi Mie 


Rearranging the above equation we get 


Me Zi Ae ey (5) 
Ig Kg Yoo} | Ye 
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Where Ic, I, and Vo, Vi indicate currents and voltages of the generator and load buses. Frc = - [Yur]! [Yu] 
are the required values. The L-index values for a given load condition are computed for all the load 
busses.The range of L-index value is 0-1. As it is closer to zero, it indicates better stability and the improved 
system security. As it approaches 1, it indicates closer to voltage collapse. Stability index L; must not be 
violated the maximum limit for any of the load buses. An L-index value away from 1 and close to zero 
indicates an improved system security. So (1-L;) indicates the margin of stability. 


2.3. +Constraints 

The two objective functions are minimized by optimally choosing the three control variables; 
Transformer tap settings, Generator excitations settings and Switchable VAR compensating settings. 

The constraints on these control variables are given as. 


iu" < ty < ty’, 1€ T 
Vim@ < Vis Vim®™, 1 E Ne (6) 
Qin < Qgi < Ogre. i Cc. Noe 


Where ti represents the tap setting of transformer connected between i-j buses, Ng is the set of generator 
buses, Vj represents the generator bus voltage of i“ bus, Q.i represents the reactive power compensation 
capacity of i™ bus and Ng is the set of load buses with reactive power support. There are two dependent 
variables, reactive power output of the generators and voltage of all load buses, which will be effected during 
optimization. So the constraints on these dependent variables need to be considered while performing 
optimization. They are given as: 


Quis Qei = Qaim, i€ Ng 
vim < Vi< ViT7@™ 1 € Ni (7) 


Qsi represents the reactive power generated by the i generator. V; is the voltage magnitude at load bus i and 
N_ is number of load buses. 

The constraints on control variables are adjusted to their limits, if they exceed, before determining 
the objective functions. The constraints on dependent variables are dealt by using penalty factor method. By 
considering this objective functions change as follows. 


ine vi-viim : (Q i O57") 
Pross = DA ge(VE — VP — 2ViV; cos(4:- 6) + DN (EO) + BUN, (ae cy) 8) 
i i gi gi 


2 ade eg 
~ N (vi-vi'™) N (Qgi- 57") 
Linax = Max(L;) + By Diy (a ae ee ih (ome onin) (9) 
By, Bz are penalty factors. V;"", Qui” can be expressed as 
VPP Ve Vet Qgi Qgi > Qi 
Vp SAV eV Oe S45 OF Oy On (10) 
V;, others Qogi, others 


3. MODIFIED BACTERIAL FORAGING ALGORITHM 
3.1. Chemo taxis insufficient of Original BFA 

In the original algorithm, search begins with population of bacteria, where each bacteria is a 
potential solution of the optimization problem. The population is converged towards optimal solution by 
following the foraging strategy of bacteria. This process consists of chemotaxis, reproduction and elimination 
and dispersion. Chemotaxis step simulates the movement of E.coli bacteria through tumbling and swimming 
via flagella. Chemotactic movement is continued until a bacteria goes in the direction of positive nutrient 
gradient i.e. increasing the fitness. It is achieved through tumbling and swimming. The chemotaxis 
movement of the bacterium can be represented as 


0G, kD = 8G kD +c 6G), 0G) = FTA (11) 
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Here 6'(j,k,1) represents the position vector of i bacterium for the j chemotaxis step, k" 
reproduction step and I" elimination dispersal step. c(i) is step size taken in the random direction specified by 
the tumble. (J) is direction angle taken by tumble at j" step. If the fitness at 6'(j + 1,k, l) is better than the 
fitness at 6'(j,k,1) then the bacterium takes another few step sizes c(i) in that direction specified by swim 
length. If the fitness at ag + 1,k,1) is not better than the fitness at 6G, k,l) then bacterium does not go for 
swim, it finds another direction through tumble. But this will affect the algorithm because many a times it 
may not find better position. It slows down the algorithm search speed and may settle at local optimum. In 
this way original chemotactic step of BFA cannot effectively stabilise the searching process. In order to 
improve the chemotaxis, different measures are to be taken to stimulate the static individual baterium to do 
extra movement. 


3.2. Differential Evolution operator 
The following differential evolution operator is selected by considering the parameter setting of 
BFA. 


vi = 04G, kD + F(09"°"G, kD) — ORG, kD + 08 G, kL) — O83, kD) (12) 


where d refers to the dimension of solution. @4,(j, k, |) represents the position vector of it" bacterium for the j" 
chemotaxis step, k'* reproduction step and 1" elimination dispersal step. v4 represents i'* bacterium during 
differential mutation . gars is the global best of the solutions, 681, @%?, 9%? are three individual bacteria 
that are randomly chosen from the whole group. F is the differential factor whose range is [0.2 — 0.9]. By 
using the information of best individual in the current population, the speed of the over all search is 
accelerated. Besides, by taking full advantage of the information of other individuals in the population, 
degradation of individual dimention is prevented and the probability of the individual trapped in to the local 
optimal is reduced. 


3.3. Incorporating Chaotic Dynamics into initial population 

Since it gives the uniform distribution function in the interval [0, 1], the tent map is chosen in 
this paper. The chaotic dynamics of the tent map is used for generating initial population. The tent map is 
defined by, 


zj** = w(1 —2|z} -0.5|) j = 12,- - -D. (13) 


where z; denotes the jth chaos variable and i denotes the chaos iteration number. Set i = 0 and generate D 
chaos variables by (@). Then i=1,2, - - - N, generate N population of initial bacteria. Then the chaos variable 
Z}. i=1,2, - - , Nis mapped in to the the search range of descision variable by the following equation. 


Xij = Xmin,j + Zi Kmail ane Lninpy = 1,2,- — —D (14) 


Where x;; is i” bacteria of j" descision variable. Xin, j 18 minimum limit on j" descision variable and Xmax, j 
is maximum limit on j" descision variable. 


3.4. Algorithmic Steps for ORPD with modified BFA 

To apply BFA algorithm, the following steps have to be followed. 
Step 1: Read the system data. Set the parameters of the BFA. 
Step 2: Choose initial population of bacteria with chaotic dynamics of tent map (13-14). 
Step 3: Elimination dispersion loop, I=1+1, k=0. 
Step 4: Reproduction loop: k=k+1, j=0. 
Step 5: Chemotaxis loop: j=j+1, Check the bacteria for the constraints. 
Step 6: Get the fitness value of objective functions (8-9) from NR load flow solution. Perform tumble by 
adding random vector to the bacteria. Calculate the fitness, if it is better than previous, perform swim for 
swim size otherwise use differential mutation operator to update positon of bacteria. If the maximum number 
of chemotactic steps (N.) is reached go to next step, otherwise go to step 5 and continue. 
Step 7: Sort the bacteria according to their fitness. Remove the worst half of the population and replace 
them with the best half. If maximum number of reproduction steps (Nre) is reached go to nextstep otherwise 
go to step 4 and continue. 
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Step 8: Eliminate the bacteria with new one with the probability of Pea i.e if a random number is greater 
than Pea. If maximum number of elimination and dispersion steps is reached go to next step otherwise go to 
step 3 and continue. 

Step 9: Print the results. 


4. RESULTS AND DISCUSSION 

The proposed modified Bacterial Foraging algorithm is applied to the ORPD problem with two 
objective functions, minimization of real power loss (Ploss) and voltage stability L-index (Vstability). The 
evolution is carried out on standard IEEE 30 bus test system. System data and initial settings are adapted 
from [17]. It consists of 30 buses, 41 branches, 6 generators, 4 tap setting transformers and 9 switchable VAR 
compensating sources. Buses 1,2,5,8,11 and 13 are generator buses. Reactive power sources are installed at 
buses 10,12, 15, 17, 20, 21, 23, 24 and 29. Branches (6-9), (6-10), (4-12) and (28-27) are equipped with 
OLTC transformers. The voltages of generator buses and load buses have been constrained within limits 
between 0.95p.u and 1.1p.u. Operating range of all OLTCs is range from 0.9 to 1.1. The range of capacitor 
banks is considered between 0 MVAr to 5 MVAr. 

Table 1 shows the simulation results for Ploss objective. The Proposed MBFA algorithm reduced 
the power loss from base value 5.812 MW to 4.5978 MW, which indicates 20% reduction from the base 
value. There is also 0.006 MW reduction of power loss in comparison with ALO algorithm which is lowest 
among other evolutionary computation algorithms presented in the table for comparison. There is also 0.11 
MW reduction of power loss in comparison with the basic BF algorithm. The proposed algorithm is also 
giving a 6% reduction of Lmax value for Pics; objective in comparison with BA which is lowest among the 
algorithms from literature. 


Table 1. Comparision of simulation results for Ploss objective 








initial BA [18] GWO[18] ABC [18] TRL pre GAFGP [20] BFA MBFA 
VG1 1.05 1.1 1.1 1.1 1.1 1.0744 1.055 1.1 1.098 
VG2 1.04 1.094 1.0938 1.0971 1.0953 1.0724 1.042 1.0956 1.0946 
VG5 1.01 1.074 1.0737 1.0866 1.0767 1.0486 1.035 1.068 1.0798 
VG8 1.01 1.076 1.0797 1.08 1.0788 1.498 1.036 1.0761 1.0817 
VGI1 1.05 1.1 1.1 1.085 1.1 1.0692 1.085 1.1 1.0965 
VG13 1.05 1.1 1.0944 1.1 1.1 1.0038 1.064 1.0953 1.1 
T6-9 1.078 0.95 0.98 1.07 1.01 1.0375 0.9536 1.0195 1.0459 
T6-10 1.069 1.03 0.97 0.95 0.99 0.9938 0.9067 0.9848 0.9052 
T4-12 1.032 0.99 1.02 1.02 1.02 0.975 0.999 1.0283 0.9759 
128-27 1.068 0.97 0.99 1.01 1 1.0438 0.9662 0.9493 0.9688 
QC10 0 5 2 5 4 0.011 0.03871 4.0169 3.5902 
QC12 0 0 5 0 2 0.033 0.04151 1.9792 4.5386 
QC15 0 5 4 2 4 0.0465 0.04812 0 3.5325 
QC17 0 3) 4 5 3 0.035 0.03735 3.0222 4.5453 
QC20 0 0 4 4 2 0.0335 0.04617 2.9253 4.8974 
QC21 0 0 0 5 4 0.018 0.04828 2.0375 1.2546 
QC23 0 0 5 4 3 0.007 0.03781 1.0387 4.4724 
QC24 0 5 a) 5 5 0.017 0.04512 4.0035 4.9146 
QC29 0 0 3 4 ) 0.0155 0.0269 2.0401 1.3196 
Ploss 4.812 4.628 4.611 4.611 4.59 5.129 5.169 4.694 4.584 
Lmax 0.1716 0.1247 0.1303 0.1326 0.1307 NR NR 0.1189 0.1176 





All the simulations are done in MATLAB R2009b software on a personal computer with 
configuration i3 processor, CPU 1.9GHz and 4GB RAM. 30 independent runs were executed and best, worst 
and mean values of optimal solutions are presented. The obtained results of proposed Modified BFA are 
compared with basic BFA and other standardard evolutionary algorithms in the literature such as Bat 
algorithm (BA), Grey wolf optimization (GWO), Artificial Bee colony (ABC), Ant Loin Optimization 
(ALO), Hybrid differential Evolution and Simulated Annealing (HDESA), Genetic Algorithm based Fuzzy 
Goal Programming (GAFGP) and Gravitational Search Optimization (GSO) algorithm. 

Table 2 the simulation results for Vstability objective. The Proposed MBFA algorithm reduced the 
power loss from base value 0.1716 to 0.1139, which indicates 33% reduction from the base value. There is 
also 2% of reduction of Lmax value from ALO algorithm which is lowest among other evolutionary 
algorithms presented in the table for comparison. The proposed algorithm also offered 3% reduction from the 
basic BFA. Figure | shows the convergence of BFA and MBFA algorithms for 100 iterations for Pos; and 
V stability Objectives. 
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Table 2. Comparision of simulation results for Vstability objective 





















































initial BA [18] GWO [18] ABC [18] ALO [18] GSO [21] BFA MBFA 
VG1 1.05 1.097 1.0965 1.0829 1.0992 il 1.0591 1.0997 
VG2 1.04 1.093 1.0807 1.073 1.0948 11 1.0569 1.0989 
VG5 1.01 1.049 1.0693 1.0759 1.0975 i 1.0409 1.0783 
VG8 1,01 1.071 1.0624 1.0744 1.0997 11 1.0893 1.0484 
VGI1 1.05 1.06 1.0977 1.1 1.0979 ile 1.0637 1.0965 
VG13 1.05 1.097 1.0927 1.0804 1.1 1.1 0.9674 1.0913 
T6-9 1.078 1.09 0.96 1.03 1.04 0.9 0.9794 0.9416 
T6-10 1.069 0.9 1.01 0.92 0.95 0.9 0.9654 1.058 
T4-12 1.032 id 0.97 0.92 0.98 0.9 0.9059 0.9884 
T28-27 1.068 0.93 0.94 0.97 0.97 1.019538 0.9325 0.9453 
QC10 0 3 2 5 = 5 3.0287 5 
QC12 0 4 1 5 3 5 3.994 2 
Qc15 0 3 1 5 3 5 3.0276 2) 
QC17 0 5 2 4 4 5 1.8953 5 
QC20 0 5 2 5 3 5 1.9395 4.9763 
QCc21 0 0 1 3 2 5 3.9801 4.9456 
QC23 0 0 4 4 1 5 4.0131 3 
QC24 0 0 4 4 2 5 4.0236 4.9768 
QC29 0 3 4 5 4 5 4.0333 4.943 
Ploss 4.812 5.0748 4.8269 4.9688 4.8693 6.660258 6.652 4.954 
Lmax 0.1716 0.1191 0.118 0.1161 0.1161 0.11607 0.1174 0.1137 
7.5 0.15 
BFA 
0.145; MBFA |7 
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Figure 1. Convergence of BFA and MBFA for Pioss and. V stability objectives 


Table 3 shows summary of the performance of both the algorithms. It clearly shows the out 
performance of the proposed algorithm for both the objectives. Generally evolutionary computation 
algorithms are random in nature they vary from one run to another run. But the standard deviation value 
shows that the proposed algorithm is more consistant which is very much desirable in practical applications. 


Table 3. Summary of BFA and MBFA performance 








BFA MBFA BFA MBFA 

Ploss best 4.694 4.584 Lmax best 0.1174 0.1137 

Ploss worst 5.138 4.706 Lmax worst 0.1258 0.1173 
Ploss mean 4.906 4.638 Lmax mean 0.1212 0.1148 
Ploss STD 0.1169 0.042 Lmax STD 0.0022 0.0011 





5. CONCLUSIONS 

Reactive power optimization with a Modified Bacterial Foraging algorithm for two objectives; 
minimization of real power loss and voltage stability L- index is proposed. The proposed algorithm is tested 
on IEEE 30 bus test system. Simulation results obtained by the proposed MBFA are compared with original 
BFA and also with other popular techinuques which are reported in the recent state of art literatures and it is 
demonstrated that there is significant improvement in both objectives in comparison with original BFA and 
also giving better results in comparison with other algorithms. The results also demonstrate that addition of 
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differential mutation and chaos dynamics improved the original BFA algoritm to a considerable degree and 
with consistency. So the proposed algorithm is suitable for Energy Control Center. 
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